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ABSTRACT 

Aims. Non-local thermodynamical equilibrium (non-LTE) line formation for singly-ionized and doubly-ionized praseodymium is 
considered through a range of effective temperatures between 7250 K and 9500 K. We evaluate the influence of departures from LTE 
on Pr abundance determinations and determine a distribution of the Pr abundance in the atmosphere of the roAp star HD24712 from 
non-LTE analysis of the Prn and Prm lines. 

Methods. A comprehensive model atom for Prii/m is presented based on the measured and the predicted energy levels, in total, 6708 
levels of Prn and Prm and the ground state of Priv. Calculations of the Prn energy levels and oscillator strengths for the transitions 
in Prn and Prm are described. 

Results. The dependence of non-LTE effects on the atmospheric parameters is discussed. At T c h > 8000 K departures from LTE lead 
to overionization of Prn and, therefore, to systematically depleted total absorption in the line and positive abundance corrections. 
At the lower temperatures, different lines of Pr n may be either weakened or amplified depending on the line strength. The non-LTE 
effects strengthen the Prm lines and lead to negative abundance corrections. Non-LTE corrections grow with effective temperature 
for the Prn lines, and, in contrast, they decline for the Prm lines. The Prn/m model atom is applied to determine the Pr abundance in 
the atmosphere of the roAp star HD 24712 from the lines of two ionization stages. In the chemically uniform atmosphere with [Pr/H] 
= 3, the departures from LTE may explain only small part (approximately 0.3 dex) of the difference between the LTE abundances 
derived from the Prn and Prm lines (^2 dex). We find that the lines of both ionization stages are described for the vertical distribution 
of the praseodymium where the Pr enriched layer with [Pr/H] > 4 exists in the outer atmosphere at log T5000 < -4. The departures 
from LTE for Prn/m are strong in the stratified atmosphere and have the opposite sign for the Prn and Prm lines. The praseodymium 
stratification analysis of roAp stars has to be performed based on non-LTE line formation. Using the revised partition function of Prn 
and experimental transition probabilities, we determine the solar non-LTE abundance of Pr as log (Pr/H) Q = -11.15 ± 0.08. 

Key words. Atomic data - Atomic processes - Line: formation - Stars: atmospheres - Stars: chemically peculiar - Stars: individual: 
HD 24712 

1 . Introduction dances derived from the Ndn and Ndm lines and not the 1.5 - 

2.0 dex observed in the these stars. Mashonkina et al. come back, 
Classical LTE analysis finds a great violation of the ionization ther ef ore , to the assumption of Ryabchikova et al. (2002) that 
equilibrium between the second and the first 10ns of the rare- the Nd anomaly observed in Equ is caused by a stratified Nd 
earth elements (REE) in rapidly oscillating chemically pecu- distribution with the accumulation of the element in the upper- 
liar (roAp) stars (Cowley & Bord [19981 Cowley et al. [20001 most atmospheric layers, above log r 50 oo = -8 according to the 
Gelbmann et al. 2000] Ryabchikova et al. [20001 Kochukhov LTE analysis It was found that ±e non . LTE effects for the Nd n 
Oj ■ [2003 Ryabchikova et_al. [200D- For the sample of 26 stars, and Nd ffl Hnes m yery stfong fa the stratified atmosphere, and 
Ryabchikova et al. <|200D show that a discrepancy between the they resuh in significant shifting the Nd enriched layer down . 
abundances derived from the lines of singly ionized and doubly ward compared t0 the location determined in the LTE analysis, 
ionized atoms of neodymium and praseodymium typically ex- The required Nd overabundance in the layer is [Nd/H] = 4 at 
ceeds 1 .5 dex in roAp stars, while it is substantially smaller if ex- j ^ < _ 3 5 for Equ and [Nd/H] = 4 5 at fo ^ < _ 4 5 
ists in non-pulsating Ap stars (see also Kato|2003] Ryabchikova for HD 247 12 
et al.[2006]for hotter Ap stars HD 170973 and HD 144897). This 

rules out errors in the oscillator strengths as a possible reason for The present paper continues to investigate the rare-earth el- 
the abundance difference observed in roAp stars. ements in stellar atmospheres based on the non-LTE line for- 
In our previous paper (Mashonkina et al. 2005 ), we inves- mation and is devoted to the praseodymium. We study the sta- 
tigate the Ndn and Ndm spectra in two roAp stars, y Equ tistical equilibrium (SE) of singly ionized and doubly ionized 
and HD 247 12, based on non-local thermodynamic equilibrium praseodymium, Pr 11 and Pr m, through a range of effective tem- 
(non-LTE) line formation and show that the non-LTE effects may peratures between 7250 K and 9500 K, evaluate the influence of 
explain only 0.5 dex of the difference between the LTE abun- departures from LTE on Pr abundance determinations in the Sun, 
A and Ap type stars, consider the non-LTE effects for Prn/in in 
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Send offprint requests to: L. Mashonkina the atmosphere with non-uniform vertical distribution of Pr, and 



finally determine empirically a stratification of Pr in the atmo- 
sphere of the roAp star HD 24712 from the non-LTE analysis of 
thePrn andPrm lines. 

The paper is organized as follows. In Sect.|2l an extensive 
model atom for Pr ii/m is introduced and theoretical calculations 
of the Pr ii atomic structure and transition probabilities for Pr n 
and Prm are presented. The programs and atmospheric models 
used in the line formation calculations are described in Sect.[3] 
Section|4]investigates the departures from LTE for Prn/m in the 
model atmospheres with homogeneous and stratified distribution 
of Pr. The non-LTE abundance corrections for the selected lines 
of Prm and Prn are given there depending on effective tem- 
perature. The solar Pr abundance is revised in Sect.[5]based on 
the improved partition function of Pr n. In Sect. [6] we determine 
the Pr abundance distribution in the atmosphere of the roAp star 
HD 247 12 from the non-LTE analysis of the Pr n and Pr m lines 
and discuss the influence of the uncertainties of atomic param- 
eters on non-LTE modelling and final results. Our conclusions 
and recommendations are given in Sect. [7] 

2. Model atom of Prn-Prm 

Model atom provides the necessary atomic input data to spec- 
ify the SE equations and the opacities/emissivities for radiative 
transfer calculations. 

2.1. Energy levels 

The lower levels in singly ionized Pr belong to the 4f 3 6s con- 
figuration with the ( 4 I) 5 I ground term. Laboratory measurements 
(Ginibre |T989alll989bl Furman et al. l200ll Ivarsson et al. [200Tb 
give 330 energy levels of Prn with an excitation energy E exc < 
5.3 eV. Most known levels belong to singlet, triplet, and quin- 
tet terms of the 4/ 3 nl (nl = 6s,5d,6p) and 4f 2 5d nl {nl = 
5d, 6s, 6p) electronic configurations. For some of the known en- 
ergy levels, only the total angular momentum and assignment to 
an electronic configuration were given. The highest known lev- 
els of Prn are separated by more than 5 eV from the ground state 
of Prm. The calculations with the Cowan (1981) code show that 
even below 5.3 eV there are many unidentified energy levels. 
With such an incomplete term system we cannot get a realistic 
statistical equilibrium of the atom. 

For the present study, we calculate the energy levels of Pr n 
using the Cowan code (Cowan 198 1 ). The measured levels of the 
odd 4/ 3 6s + 4/ 3 5ii + 4f 2 5d6p configurations have been fitted 
taking into account the interactions with the 4f 2 6s6p, 4-f5d 3 , 
4f5d 2 6s, and 4f5d6s 2 configurations. The 4f 2 5d 2 , 4f 2 5d6s, 
4f 2 6s 2 , 4/ 3 6p, 4/ 4 , 4f 2 6p 2 , and 4f5d 2 6p configurations are 
included in the fitted matrix of the even system. The fitting re- 
sults in a standard deviation of the calculated from experimen- 
tal levels of 86 cirT 1 for 129 odd and 72 cirT 1 for 201 even 
levels. To get highly excited levels the energy structures of the 
even 4f 3 (lp - 9p) and 4/ 3 (5/ - 6f) configurations and the odd 
4/ 3 (7s - 9s), 4/ 3 (6<i - Sd), 4/ 3 5g, and 4f 2 5d5f configurations 
are calculated. In all unknown configurations, the average en- 
ergies where Brewer (119711) predictions are absent and the en- 
ergy parameters of the 4/ 2 and 4/ 3 cores are scaled similarly to 
the known configurations, the other parameters being scaled by 
a factor 0.75 with respect to the corresponding Hartree - Fock 
values. As a result, a nearly complete set of the levels below 
70000 cirT 1 (8.7 eV) has been obtained consisting of 3938 ener- 
gies. All the known levels and the predicted levels with E exc < 
10.04 eV, in total 6539 levels of Prn, are used to construct the 
model atom. They are shown in Figs.Q] and [2] (Online only) for 



Table 1. Partition functions for Prn and Prm. 



Temperature 


Prn 


Prm 


(K) 


Ginibre (1989) 


this study 


this study 


3000 


55 


58 


22 


4000 


91 


101 


28 


5000 


138 


166 


34 


6000 


192 


251 


43 


7000 


252 


356 


54 


8000 


315 


478 


66 


9000 


380 


617 


81 


10000 


446 


770 


97 


11000 


511 


937 


114 


12000 


576 


1117 


132 


13000 


639 


1309 


151 


14000 


701 


1514 


170 


15000 


762 


1732 


190 


16000 


820 


1962 


211 


17000 


878 


2205 


231 


18000 


933 


2460 




19000 


986 


2727 


273 


20000 


1038 


3007 


293 



the even and odd levels, correspondingly. The calculated high 
excitation levels provide the close collisional coupling of Pr n to 
the continuum electron reservoir. 

It is worth noting that the calculations lead to a significantly 
larger partition function of Pr n compared to that based on the 
laboratory levels only. This is illustrated in TableQ] where the 
partition function of Pr n is calculated by its definition for a tem- 
perature range 3000 K - 20000 K using the measured energy lev- 
els from Ginibre ( 1989a 1989b) and using the levels predicted in 
this study. In Sect. [5] we calculate the effect of the revised parti- 
tion function on the Pr abundance determined from the solar Pr n 
lines. 

For Prm, the laboratory measurements resulted in 593 en- 
ergy levels (Martin et al. 119781 Palmeri et al. 120001) with an ex- 
citation energy up to 17.5 eV. In the range of stellar parameters 
we are concerned with, there is no need to include the highly ex- 
cited levels of Pr m in the model atom. They play a minor role in 
population and depopulation of Pr m, because the next ionization 
stage Pr iv represents a negligible fraction of Pr abundance. We 
use, therefore, the levels of the odd 4/ 3 and the even 4f 2 5d and 
4f 2 6s electronic configurations with E exc < 6.7 eV. The term 
structure is shown in Fig. [3] (Online only). The contribution of 
the Prm energy levels omitted in the final model atom to the 
Prm partition function is less than 0.01% at T e = 7250 K. 

Levels of the same parity with small energy differences were 
combined into a single level. The final model atom includes 294 
combined levels of Prn, 54 combined levels of Prm, and the 
ground state of Pr iv. 

2.2. Radiative data 

The 15788 and 392 radiative bound-bound transitions in Prn 
and Pr in, respectively, are included in SE calculations. Oscillator 
strengths fij based on laboratory measurements of Lage & 
Whaling (1976 ) or estimated from the observed line intensities 
of Meggers et al. ( 119761 ) are available only for 448 transitions in 
Prn. They come from the Kurucz & Bell ( 1995) linelist and are 
accessible via the Vienna Atomic Line Data Base (vald, Kupka 
etal. [19991 . 

For the majority of transitions in Prn and all transitions 
in Prm, we rely, therefore, on oscillator strengths computed in 



the present study. Calculations for Pru are based on the wave 
functions obtained in the fittings of the energy levels. All the 
Hartree-Fock transition integrals are scaled by a factor 0.85. The 
calculated lifetimes have been compared with the recent accurate 
measurements of Scholl et al. (2002) and Biemont et al. (2003). 
A quite good agreement is found for the levels with leading 
contribution from the levels of the 4/ 3 6p configuration: the ratio 
of the measured to calculated lifetime equals 1 . 1 with a standard 
deviation 0.3. For the levels with highly mixed wave functions, 
the corresponding quantity appears to be 0.7 with the standard 
deviation 0.5. Details of these calculations can be found on 
http: //das 1(91 .isan.troitsk.ru/files/SPECTRA/Pr_II 
and will be presented in a forthcoming paper. 

For Prm, our approach follows that of Palmeri et al. (2000). 
Fitting of the energy levels is performed taking into account 
the interactions in the odd 4/ 3 + 4f 2 6p + 4f5d 2 + 4f5d6s 
+ 4f6s 2 + 4/ 2 5/ + 5pHf 4 and the even 4/ 2 5«i + 4f 2 6d 
+ 4/ 2 6s + 5pHf5d + 5pHf6s + 4f5d6p + 4f6s6p + 
4f5d5f complexes. Ab initio transition integrals are taken. 
For 7 levels with a more than 25% contribution of the 4f 2 6p 
configuration, the average ratio of the measured (Biemont et al. 
1200 11 1 to calculated lifetimes is obtained to be 1.17 + 0.25. 
An exception is the J = 9/2 level at 62535.6 cirT 1 with 
main contribution from the 4f5d 2 configuration and about 
11% admixture of the 4f 2 bp configuration. For this level, 
the lifetime ratio appears to be 3.1. Details of the Prm cal- 
culations relevant to this article can be found on the website 
http : //das 1(8 1 . isan . troitsk . ru/f iles/SPECTRA/Pr_III 

The photoionization cross-sections o" p h are computed using 
the hydrogen approximation because no accurate data is avail- 
able for the Pru and Prm levels. We assume that the pho- 
toionization from any Pru level ends in the ground state of 
Prm, 4f 34 Igj T To take into account the photoionization to the 

remaining levels of the 4 _f 3 4 I° term we multiply <r p h by the 
ratio g(4f 34 l°)/g(4f 34 r 9/2 ) = 5. The photoionization from the 

4f 2 5dnl levels ends in the excited Pr m 4f 2 5d levels, however, 
we neglect this for the following reason. Replacing the real end- 
ing state on the Pr m ground state makes each 4f 2 5dnl level to 
be easier ionized. Its ionization energy is reduced by more than 
2 eV, and the photoionization rate is overestimated, in particu- 
lar, for the levels with a threshold in the ultraviolet (UV). At the 
same time, the photoionization rates of the 4f 2 5dnl levels are 
underestimated due to ignoring multiple channels for their ion- 
ization. We cannot estimate net effect due to absence of accurate 
atomic data. 

2.3. Collisional data 

The calculations of electron impact excitation and ionization 
rates rely on theoretical approximations. We use the formula of 
van Regemorter (1962 ) for the allowed transitions and assume 
that the effective collision strength T = 1 for the forbidden tran- 
sitions. Electron impact ionization cross-sections are computed 
according to Drawin d 1 96 It . 

The effects of the uncertainties of the used photoionization 
cross-sections and collisional rates on the final results are de- 
scribed in Sect. 16.31 

3. Programs and model atmospheres 

The radiative transfer and the statistical equilibrium equations 
are solved with a revised version of the DETAIL program (Butler 
& Giddings 1985 ) using the accelerated lambda iteration fol- 



lowing the extremely efficient method described by Rybicki & 
Hummer (19911 [1992) . Background opacities include the im- 
portant bound-free and free-free transitions of hydrogen, he- 
lium, and the most abundant metals, the Rayleigh scattering, 
the Thomson scattering, the hydrogen lines, the quasi-molecular 
Lyman a satellites due to H-H and H-H + collisions (Allard 
et al.1998 ), and line opacity calculated with the line lists made 
available by Kurucz & Bell (1995 ). The background opacities 
are sampled on a random grid of 4500 frequencies, to which 
are added the frequencies of the line profiles. The final non-LTE 
line formation program samples roughly 150 000 wavelengths 
between 500A and 80 000A. 

The obtained non-LTE and LTE level populations are used to 
compute the emergent flux, line profiles, and equivalent widths 
with the code LINEC. The investigated lines and transitions of 
Pr n and Pr m are listed in Table|2] The wavelengths of Pr n lines 
are taken from Ginibre ( 1990 ). Ritz wavelengths are used for the 
Prm lines. 

Hyperfine-structure (HFS) affecting the Pr n lines is explic- 
itly calculated with the HFS constants given by Ginibre ( 1989a), 
if available. No HFS data exist in the literature for the Pr m lines. 
Ivarsson et al. (2001 , ILW) could not resolve HFS in any of the 
accurately measured Prm lines, which means that HFS is not 
large, although some lines show asymmetric profiles with half- 
widths exceeded the half-width of the resolved HFS components 
in Pr n. The three from eight common Pr m lines are practically 
HFS unaffected according to the ILW laboratory analysis. Half- 
widths of other three lines are three to four times larger com- 
pared to those for the HFS unaffected lines, and they are com- 
parable with the thermal half- widths of the Prm lines in the at- 
mosphere of HD 24712. The thermal half-width of the Pr lines 
corresponds to approximately 1 kms _I . Our simulation of HFS 
affecting the Prm 5998 and 6053 lines with the strongest HFS 
broadening according to the measurements of Ivarsson et al. 
(2001 ) leads to an increase of the theoretical equivalent widths, 
by 16 mA and 17 mA, respectively, that transforms to -0.20 dex 
and -0.24 dex decrease of the element abundance derived from 
these lines. Further analysis is performed not accounting for HFS 
for the Prm lines. We show in Sect. 16.21 that two abovemen- 
tioned lines give the results consistent with those for the Prm 
lines with weak HFS broadening. 

We use homogeneous blanketed model atmospheres. The 
small grid of models with r e jf ranging between 7500 K and 
9500 K with a step of 500 K, log g = 4, and the solar chem- 
ical composition has been calculated with the MAFAGS code 
(Fuhrmann et al. l 19971 ) that treats line-blanketing using the opac- 
ity distribution functions (ODF). The original ODF tables of 
Kurucz ( 1994a) were scaled by -0.16 dex to put the iron opacity 
calculated by Kurucz with logep e = 7.67 into correspondence 
with a value logSFe =7.51 which we believe to be the best rep- 
resentation of the solar mixture. We refer to abundances on the 
usual scale where log su = 12. 

For HD 24712 (T eff = 72 50 K log g = 4.3, [M/H] = ac- 
cording to Ryabchikova et al. 11997) , the model atmosphere has 
been computed by Frank Grupp with the MAFAGS -OS code 
(Grupp 2004). It is based on up-to-date continuous opacities and 
includes the effects of line-blanketing through opacity sampling. 

HD24712 possesses a magnetic field with a mean magnetic 
field modulus (B) changing from 2.5 kG at the magnetic mini- 
mum to ~3. 1 —3.3 kG at the magnetic maximum (Ryabchikova 
et al. 12007b) . We ignore the influence of the magnetic field on 
atmospheric structure, based on the results of Kochukhov et al. 
(2005 ) who have shown that the difference in temperature and 



gas pressure distributions between magnetic and non-magnetic 
model atmospheres with T e ff = 8000 K does not exceed 30 K and 
6%, respectively, for field strengths up to 5 kG. However, split- 
ting of the spectral lines in the magnetic field is taken into ac- 
count in LTE abundance analysis. Magnetic spectrum synthesis 
is performed with the help of SYNTHMAG code (Kochukhov 
120071 . 

There are observational evidences for non-uniform element 
distribution in the atmosphere of HD 24712. The neodymium is 
strongly enhanced in the uppermost atmospheric layers accord- 
ing to Mashon kina e t al. (120051 . Ryabchikova et al. (120081 1 and 
Ryabchikova (2008) show that Ca, Si, Cr, Fe, Sr, and Ba are con- 
centrated in deep atmospheric layers. One believes that radia- 
tively driven diffusion is one of the main processes responsible 
for these inhomogeneities. The self-consistent diffusion models 
predict the vertical distributions of Mg, Si, Ca, Ti, Fe which 
qualitatively reproduce the corresponding element stratifications 
found empirically for some Ap stars (Alecian & Stift [2007 
LeBlanc & Monin 2004J. But no theoretical predictions are 
available for the REE due to an incompleteness of atomic data 
on energy levels and transition probabilities. Stratified distribu- 
tion of chemical elements in the atmosphere can influence the at- 
mospheric structure. Modelling chemically non-uniform stellar 
atmospheres based on the abundance gradients determined em- 
pirically is in progress (Shulyak, 2008, private communication). 
Preliminary results of the iterative procedure for HD 24714 pre- 
dict the slight change in the model parameters (T e s = 7250 K 
and log g = 4. 1) and lead to only minor changes in our results on 
Pr analysis presented in Sect.|6] 

4. Departures from LTE for Prii-Prin 

In line formation layers of the atmospheres with r e fj between 
7250 K and 8000 K, number density of Pr n is larger compared 
to that for Prm (see Fig.|U. However, Prn drops rapidly in the 
hotter atmospheres. Non-LTE calculations show that the main 
non-LTE effect for Prn is overionization caused by a super- 
thermal radiation of non-local origin near the thresholds of the 
4f6p levels with £ exc = 3 eV - 4 eV (A th[ = 1600A to 1850A). 
Photoionization is able to drain the populations of these levels 
in the model with the lowest temperature, T e ff = 7250 K, and 
the effect is strengthened with T e ff increasing. The population 
loss is redistributed over many levels, producing overall deple- 
tion of the first ionization stage in line formation layers. In con- 
trast, photoionization of the Pr in levels is inefficient due to very 
low stellar fluxes in the far ultraviolet (A < 900A), where the 
ionization edges of the Prm ground state and the low excita- 
tion levels are located. At T e g < 8000 K, overionization of Pr n 
leads to overpopulation of Pr m. At the higher temperatures, Pr in 
represents the majority of the element and preserves the thermo- 
dynamic equilibrium (TE) total number density. 

We find that the mechanisms driving departures from LTE 
for Pr ii/iii are similar in the atmospheres with uniform and strat- 
ified distribution of praseodymium, however, the magnitude of 
the effect is different due to a different location of the line for- 
mation layers. 

4.1. The atmospheres with uniform distribution of 
praseodymium 

In the observed spectrum of HD 24712, 14 lines of Prm are de- 
tected with a measured equivalent width, Wa, of 17 mA and 
larger. In order to predict the equivalent widths of the Prm 
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Fig. 4. LTE and non-LTE total number densities of Pr n and Pr m 
(top panel) and departure coefficients for selected levels of Pr n 
and Prm (bottom panel) in the model atmosphere 7250/4.3/0. 
Successive numbers of the levels in our model atom are quoted. 
Everywhere in the atmosphere [Pr/H] = 3. 



lines at a detectable level in the model, which represents the 
atmosphere of HD 24712, we perform calculations with the 
praseodymium abundance [Pr/H] = 3. This value also charac- 
terizes the mean Pr abundance derived from all investigated Pr n 
and Prm lines in HD 24712. FigureH] shows the departure coef- 
ficients, bj = n^ LTE /ti™ of the selected levels of Prn and Prm 

as a function of continuum optical depth T5000 at A = 5000A 
in the model atmosphere with r eff /log g/[M/H] = 7250/4.3/0. 
Here, «^ LTE and n™ are the statistical equilibrium and TE 
(Saha-Boltzmann) number densities, respectively. The medium 
becomes optically thin for the ionizing radiation below 1850A 
far inside the atmosphere, at log T5000 around 0. As a result, Pr m 
is overpopulated in line formation layers, at log T5000 < 0. The 
Prn levels with E exc < 6 eV are strongly coupled to the Prn 
ground state and to each other inside logrsooo - —1.7 where 
their departure coefficients are only slightly below 1. The out- 
side layers become transparent for the radiation of many Prn 
lines arising between low-excitation terms (E exc < 2 eV) and 
intermediate-excitation terms (E exc = 2.7 - 4 eV). The photon 
loss in these lines amplifies the underpopulation of the upper 
levels (level numbers from 110 to 146 in Fig. 01 caused by en- 
hanced photoionization. 

The non-LTE effect on the line strength is determined by the 
departures from LTE for the lower and upper levels of the tran- 
sition in the line formation layers. For the model 7250/4.3/0, the 
non-LTE effects are very small for the Prn lines. The obtained 



Table 2. Investigated transitions in Pr n and Pr m. 



A [A] 


^''exc \P * ) 


Transition 






Pru 










4222.93 1 * 


0.05 


4/ 3 5rf 5 L>4/ 3 6p 5 K 6 


2- 


118 2 


4449.83 1 * 


0.20 


4f 3 5d 5 l° 6 -4f 3 6p 5 K 6 


3 


118 


5002.44 


0.80 


4f5d 5 iq-4f6p 3 K 6 


12 


- 129 


5110.76' 


1.15 


4f5d 5 L° m -4f6 P s K 9 


26 


- 143 


5129.54* 


0.65 


4/ 3 5tf 5 K°-4/ 3 6p 5 I 5 


10 


- 121 


5135.14* 


0.95 


4f 3 5d 5 K; -4/ 3 6/7 5 K 8 


16 


- 134 


5259.73 1 * 


0.63 


4/ 3 5tf s L°-4/ 3 6p s K 6 


9 


118 


5292.62* 


0.65 


4/ 3 5rf s K°-4/ 3 6p 5 K 6 


10 


- 118 


5322.77* 


0.48 


4/ 3 5rf 5 L°-4/ 3 6p 5 K 5 


6 


110 


5681.88 


1.16 


4/ 3 5d 5 H^ -4/ 2 5rf 23 G 5 


26 


- 134 


6017.80* 


1.11 


4/ 3 5</ 5 G°-4/ 3 6p 5 H 3 


23 


- 126 


6165.94* 


0.92 


4/ 3 5rf s I°-4/ 3 6 P 5 I 4 
4/ 3 5rf 3 L°-4/ 3 6p 3 K g 


15 


- 116 


6656.83* 


1.82 


60 


- 146 


Prm 








4910.82 


0.17 


4/ 34 In /2 -4/ 2 5rf 4 H n/2 


296 


-322 


4929.12 


0.36 


4/ 34 I, 3/2 -4/ 2 5rf 4 H 13/2 


297 


-325 


5284.69 


0.17 


4/ 34 I 11/2 -4/ 2 5rf 4 H 9/2 


296 


-319 


5299.99 


0.36 


4/ 34 F, 3/2 -4/ 2 5rf 4 H u/2 


297 


-322 


5844.41 


1.24 


4/ 32 H^ /2 -4/ 2 5rf 2 G 7/2 


300 


-332 


5998.97 


0.17 


4/ 34 In /2 -4/ 2 5rf 4 G 9/2 


296 


-315 


6053.00 


0.00 


4/ 34 i; /2 -4/ 2 5rf 4 G 7/2 


295 


-311 


6090.01 


0.36 


4/ 34 I 13 /2-4/ 2 5rf 4 H u/2 


297 


-318 


6160.23 


0.17 


4/ 34 I 11/2 -4/ 2 5rf 4 H 9/2 


296 


-314 


6195.62 


0.00 


4/ 34 I° /2 -4/ 2 5d 4 H 7/2 
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Fig. 5. Non-LTE abundance corrections for the Pru and Prm 
lines depending on T e s- For all models, log g = 4. The calcu- 
lations were made with the Pr abundance [Pr/H] = 3. 
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overpopulation of the Pr m levels leads to the strengthening of 
the Prm lines compared with the LTE case. The theoretical non- 
LTE and LTE equivalent widths of the selected lines and the 
non-LTE abundance corrections Anlte = log £nlte - log slte are 
presented in Table|4](columns 6 - 8). Oscillator strengths used in 
these calculations are taken from Kurucz & Bell ( 1 19951 1 for the 
Pr ii lines and are computed in this paper for the lines of Pr m. 

As expected, the departures from LTE for the Pr n lines grow 
with r e ff. The larger T e ff, the stronger ultraviolet radiation is, 
thus resulting in amplified overionization of Pr n. Figure[5]shows 
the calculated non-LTE abundance corrections for some lines of 
Pru and Prm depending on T e g. Everywhere, [Pr/H] = 3 is 
adopted. For the Pru lines, Anlte are positive. They are small 
at r eff = 7250 K and 7500 K with A NLTE < 0.1 dex and grow 
to approximately 0.4 dex at T s s - 9500 K. In the hottest model, 
Pr ii /15002 and /16656 have W,\ < 1 mA, and their Anlte are not 
shown. In contrast, the departures from LTE for Pr m are weak- 
ened toward higher T e ff because doubly-ionized praseodymium 
tends to represent the majority of the element and to preserve the 
TE total number density. The non-LTE abundance corrections 
are negative for the Prm lines with Anlte - (-0.3) - (-0.4) dex 
at r eff = 7500 K and A NLT e - -0.1 dex at 7/ eff = 9500 K. 

In this paper, we present also the revised non-LTE abundance 
corrections for the Nd ii and Nd m lines (Fig.|6]l. Opacity package 
of the DETAIL code was recently updated by the inclusion of the 
quasi-molecular Lyman a satellites following Castelli & Kurucz 
(2001) implementation of the Allard et al. ( 1998) theory and by 



Fig. 6. Non-LTE abundance corrections for the Nd n and Nd m 
lines depending on T e g. For all models, log g = 4. The calcula- 
tions were made with the Nd abundance [Nd/H] — 2.5. 



the use of the extended line list based on not only measured but 
also predicted lines from Kurucz & Bell ( 11995I I. The effect of the 
increased UV opacity below 1700 A is seen only for the lines of 
Nd ii at r e ff > 8500 K. At the lower temperatures, the revised 
Anlte of the Ndn lines agree within 0.01 dex to 0.03 dex with 
that obtained by Mashonkina et al. (1200 5). For the Ndm lines, 
the maximal difference between the revised and 2005 's Anlte 
constitutes 0.05 dex at T e ff = 8500 K and is smaller at the lower 
and higher temperatures. Thus, the Nd abundance distributions 
in the atmospheres of HD 24712 (T e ff = 7250 K) and y Equ {T e g 
= 7700 K) found in our earlier paper do not need to be revised. 
For the Nd n lines at T e ff > 8500 K, we obtain now the smaller 
departures from LTE due to decreased ionizing radiation. The 
difference in Anlte equals 0.1 dex to 0.2 dex for different lines 
at different temperatures. 

4.2. The atmospheres with stratified distribution of 
praseodymium 

In this section, the departures from LTE for Pr ii/m are investi- 
gated in the model 7250/4.3/0 with the stratified Pr abundance 
distribution determined in Sect. 16.21 and shown by continuous 
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Fig. 7. Top panel: the [Pr/H] ratio distributions found in the 
atmosphere of HD24712 (continuous line) and from the test 
calculations with 2 dex smaller photoionization cross-sections 
(dashed line) and with a variable effective collision strength 
for forbidden transitions (dotted line). Bottom panel: the depar- 
ture coefficients log b for the selected levels of Pr n and Pr m in 
the model 7250/4.3/0 representing the atmosphere of HD 24712 
from calculations with the praseodymium distribution shown in 
the top panel by continuous curve. See text for more details. 



curve in the top panel of Fig. [7] The departure coefficients for the 
selected levels of Pru and Prm are shown in the bottom panel 
of the same figure. The theoretical non-LTE and LTE equivalent 
widths together with non-LTE abundance corrections are given 
in Table [4] (columns 9 - 11). 

Non-LTE calculations show the depletion of Pru and en- 
hanced number density of Prm in the line formation layers 
similar to the case of homogeneous Pr abundance distribution. 
However, three important distinctions in the behavior of depar- 
ture coefficients can be seen. 

- The magnitude of the depletion of Pru is much larger in 
the stratified atmosphere, in particular, above log T5000 = - 4 
where the lines are formed. 

- In the stratified atmosphere, for every pair of the Pr 11 levels 
with i < j, an inequality bi < bj is valid. 

- Above log T5000 = — 4 in the stratified atmosphere, the excited 
levels of Pr m are decoupled to the ground state to the more 
extent compared to the case of the homogeneous Pr abun- 
dance distribution. 

In the stratified atmosphere, the Pr lines are formed in the 
uppermost atmospheric layers, above log T5000 = - 4, where the 



departures from LTE are large. Non-LTE leads to weakening the 
Pr 11 lines due to the decreased number of absorbers (bi < 1) and 
due to the line source function Stj - b j/bjB v (T e ) exceeding the 
Planck function (bj < bj). The non-LTE abundance corrections 
for various Pru lines are positive and may reach +1.2 dex. The 
Prm lines are strengthened compared with the LTE case, and 
Anlte may reach -0.7 dex. Since the non-LTE effects have the 
opposite sign for the Pr 11 and Pr in lines, they are important for 
the comparison of Pr abundances derived from these lines. 

5. Solar abundance of praseodymium 

We first apply the non-LTE calculations to the Pr abundance 
analysis for the solar atmosphere. The earlier determinations 
were performed by Biemont et al. (1 19791 1 based on the nine lines 
of Pru in the 3990 - 5330 A spectral region and by Ivarsson 
et al. (120031 1 based on the three lines of Pru. The obtained 
mean abundances of the praseodymium are surprisingly differ- 
ent, log (Pr/H) = -11.29+0.08 and log (Pr/H) = -11.60+0.1, 
respectively. The most recent study is based on the modern accu- 
rate laboratory measurements for oscillator strengths of Ivarsson 
et al. (1200 11 1. Ivarsson et al. (2003) also test the influence of 3D 
effects on the determination of the Pr abundance and draw a pre- 
liminary conclusion that the differences in the equivalent widths 
of the Pru lines between the 3D and ID cases are insignificant. 

Each of the lines used by our predecessors was checked for 
blending using the NSO solar flux spectrum (Kurucz et al. 1984) 
observed with the spectral resolving power R st 340000 at wave- 
lengths between 4000A and 4700A and with R = 520000 at 
longer wavelengths. We find four lines suitable for precise spec- 
tral fitting. They are listed in Table [3] Spectral region around 
each investigated line was synthesised with the SIU code (Reetz 
119911 1. For given atomic level, SIU computes the non-LTE pop- 
ulation as the production of the LTE occupation number and the 
corresponding departure coefficient. In LTE calculations, the re- 
vised Pr 11 partition function is applied. Calculations of the Pr 11 
lines were made with the three different sets of transition proba- 
bilities taken from Kurucz & Bell ( 1995), Ivarsson et al. ( 120011 1. 
and Li et al. (120071 1. Atomic parameters for other atomic lines in 
the synthesised regions are taken from the latest release of VALD 
(Kupka et al. 1999J. For molecular lines, we apply the data 
compiled by Kurucz dl994bl >. We use the theoretical ATLAS9 
model atmosphere of the Sun (5777/4.44/0, Heiter et al. 120021 
with convection treated according to Canuto et al. ( 1 19961 1 and 
the semi-empirical model of Holweger & Miiller (119741 HM). A 
microturbulence velocity V m i c = 0.9kms~' is adopted. Our syn- 
thetic flux profiles are convolved with a profile that combines a 
rotational broadening of 1.8 kms -1 , broadening by macroturbu- 
lence with a radial-tangential profile of V mac - 3.6 km s _1 , and in- 
strumental broadening with a Gaussian profile corresponding to 
the spectral resolution. A macroturbulence velocity was allowed 
to vary within 0.2kms~' to achieve the best fit to the observed 
line shape. 

The best fits to the observed solar line profiles are shown 
in Fig. |8] The results from calculations with the ATLAS9 solar 
model atmosphere are presented in Table[3] Below we give brief 
notes on the individual lines. 

Pr 11 4222.93 A. The continuum flux at 4222.93 A is influ- 
enced by the far wing of the Ca 1 resonance line at 4226.7 A. 
The Cai line is treated using the measured oscillator strength 
from Smith & Gallagher (1966) and the van der Waals broad- 
ening parameters based on the advanced perturbation theory of 
Anstee & O'Mara JI995I The Pru 4222.93A line is partially 



ovelapped with the blue wing of the spectral feature produced 
by two CH lines at 4223.091 A and 4223. 1 13A. A half-width of 
this spectral feature is only slightly affected by the Pr n line and 
is fitted reasonably well if we reduce oscillator strengths for both 
CH lines by 0.3 dex compared to those given by Kurucz ( 1994b). 
Small changes in wavelengths (no more than 0.014 A) and oscil- 
lator strengths (no more than 0.27 dex) were introduced to fit the 
Cr 1 4222.740 A and Fe 1 4223.237 A lines just to make better fit. 

Pru 5259.7 A. This line is well isolated and very well fitted 
by spectrum synthesis. The Ni i 5259.466 A and Ti i 5259.973A 
lines shown in Fig.[8]do not affect the Pr abundance determina- 
tion. For the Ni i line, changes in wavelength by 0.01 1 A and in 
oscillator strength by -0.4 dex were made to fit the line profile, 
while a 0.005 A wavelength change is required for Ti i line. 

Pru 5322.8 A is well isolated, too, but the continuum flux is 
influenced by the far wing of the Fe i 5324. 179 A line. With the 
best atomic parameters of the Fei line, log gf = -0.103 (Bard 
et al. [T99~n i and logy 6 /N H = -7.035 at T e = 10000 K (inc reased 
by 0.2 dex compared with the Anstee & O' Mara's 1 1 9951 value), 
the predicted synthetic spectrum is still 0.2% higher than the ob- 
served spectrum around 5323 A. The local continuum level was, 
therefore, adjusted to fit the blue wing of Fe i 5324 A. A varia- 
tion of 0.2% in the continuum flux produces a 0.05 dex change 
in the Pr abundance derived from Pru 5322.8 A. 

Pr ii 4449.83 A is located between the Dyn 4449.70 A and 
Si 1 4449.90 A lines. The blue wing and the core of the Dy n line 
are slightly affected by the adjacent lines and are well fitted using 
a center line wavelength of 4449.707 A, accurate experimental 
oscillator strength, log gf - -1.03 (Wickliffe et al. 2000}, and 
the solar Dy abundance logepy,© = 1-14 (Lodders 2003). The 
Si i line is predicted, and its oscillator strength has to be reduced 
down to log g f — -3.33 to fit the line core that is only slightly af- 
fected by the HFS components of the Pr n line. In Fig.[HJ we show 
the best fit of the 4449 A blend and also the synthetic spectra 
computed with +0.05 dex change in adopted Pr abundance. The 
uncertainty of the element abundance derived from Pr n 4449 is, 
obviously, does not exceed 0.05 dex. 

In the solar atmosphere, Pru represents the majority of the 
element, and its ground state and the low excitation levels which 
are the lower levels of the investigated transitions preserve the 
TE number density. Our non-LTE calculations show that pho- 
ton pumping produces enhanced excitation of the upper levels 
resulting in weakening the lines of interest compared to the LTE 
case. We find that the statistical equilibrium of Pru in the solar 
atmosphere is sensitive to a variation of collision excitation rates 
which include, for the Sun, interactions with not only electrons 
but also neutral hydrogen atoms. For hydrogenic collisions, we 
use the formula of Steenbock & Holweger (1984) for allowed 
transitions and follow Takeda (1994) for forbidden transitions. 
Both theoretical approximations provide only an order of mag- 
nitude estimate. Empirical constraining the efficiency of hydro- 
genic collisions in the SE of various atoms that is represented 
by a scaling factor Sh applied to the mentioned formula gives 
a variety of estimates between Sh = 0.002 and Sh = 1 (for re- 
view, see Mashonkina 2008). When pure electronic collisions 
are taken into account (Sh — 0), the non-LTE abundance correc- 
tions constitute +0.03 dex to +0.08 dex for different lines. The 
non-LTE effects become negligible (Anlte ^ 0.01 dex) when 
both types of collisions are included and Sh = 1- We adopt an 
intermediate value, Sh = 0.1. The corresponding non-LTE cor- 
rections are shown in the last column of Table [3] 



The average LTE and non-LTE praseodymium abundances 
from calculations with the ATLAS9 solar model atmosphere and 
their standard deviations are presented in two bottom strings 
of Table [3] The Pr abundance determined with the HM model 
is larger, by 0.05 dex, on average. It is worth noting that the 
statistical error of our abundance determination estimated from 
the line-to-line scatter corresponds, in fact, to the quoted er- 
rors of logg/ measurements which varies between 0.03 dex 
and 0.07 dex for the lines in Table [3] The three solutions found 
with log gf from different sources are consistent within the er- 
ror bars. However, the preference should be given to the ILW 
and LCH abundances based on the modern experimental tran- 
sition probabilities. The average ILW and LCH abundances 
agree within 0.01 dex. For each set of atomic parameters and 
for both solar model atmospheres, the derived Pr abundance 
is 0.09 dex - 0.23 dex larger than that obtained by Biemont 
et al. i 19791 1 and 0.40 dex - 0.54 dex larger compared to that 
of Ivarsson et al. (2003). Different studies give the meteoritic CI 
Chondrites abundance log(Pr/H) met = -11.20+0.04 (Grevesse 
et al. [19961 Palme & Jones ISjU^-L 1.22+0.03 (Lodders HH 
or -11.25±0.03 (Asplund et al. [2D05b . Our best solar Pr abun- 
dance for the ATLAS9/HM model, log (Pr/H) = -11.15 + 0.08 
/ — 1 1 .10+0.08, exceeds the meteoritic values by 0.05 dex to 
0.10 dex. 

In solar abundance determinations, we use the partition func- 
tion calculated by its definition and based on an extensive set of 
the predicted energy levels of Prn, t/ ca / c (Prn). As can be seen 
from TableQ] £/ ca / c (Prn) is about 30% larger than that for the ex- 
perimental energy levels only, U exp (Pr n), at temperatures of the 
line formation layers in the solar atmosphere. 

With C/ fvp (Pni), the inferred solar Pr abundance is 
log(Pr/H) = -11. 22 and -11.17, for the ATLAS9 and the 
Holweger & Miiller d 19741 ) models, respectively. Thus, accurate 
solar Pr abundance depends on the completeness of known Pr n 
levels and on the accuracy of their energies. It seems that the 
current theoretical energy calculations give us the upper limit 
for the solar Pr abundance. 

6. Praseodymium in the roAp star HD24712 

We study the Pr abundance in the atmosphere of the roAp star 
HD24712 using its average spectrum obtained during spectro- 
scopic monitoring on November 11/12, 2004 with the UVES 
spectrograph at the 8.2-m telescope, UT2 (Kueyen), of the VLT 
on Paranal, Chile (programme 274.D-5011). Observations and 
data reduction are described by Ryabchikova et al. (2007b). The 
spectrum has been obtained near the magnetic maximum, that 
coincides roughly with the REE spectral line intensity maxi- 
mum. The surface magnetic field measured at the magnetic max- 
imum is (B) = 3.1 - 3.3 kG (Ryabchikova et al. l2007bb . 

6.1. The problem of the Pr abundance in HD 24712 

We start from element abundance determination using the lines 
of Pru and Prm listed in Tableland assuming the uniform dis- 
tribution of Pr in the atmosphere of HD 24712. The choice of 
the spectral lines was determined mainly by the possibility of 
pulsation measurements in the HD24712 spectrum. Magnetic 
spectrum synthesis is performed with the help of SYNTHMAG 
code under the LTE assumption using the model atmosphere 
7250/4.3/0. We apply most recent experimental transition prob- 
abilities by Ivarsson et al. (120011 1 and Li et al. (120071) if available 
and oscillator strengths from Kurucz & Bell ( 1 19951 ) for the re- 
maining Prn lines. The derived LTE abundance is log (Pr/H) = 



Table 3. Solar praseodymium abundance from calculations with the ATLAS9 solar model atmosphere. 
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4222.93 


0.13 


0.27 


0.24 


-11.14 


-11.28 


-11.25 


+0.03 


4449.83 


-0.32 


-0.17 


-0.26 


-11.07 


-11.22 
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0.08 


0.11 
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-11.10 
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-0.32 
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-11.04 
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+0.02 






LTE-average 




-11.09 + 0.04 


-11.17 ±0.10 


-11.18 ±0.08 








non-LTE-average 




-11.06 ±0.04 


-11.14 ± 0.10 


-11.15 ±0.08 





1 Kurucz& Bell (1995 , KB) 2 Ivarsson et al. (2001 . ILW) 3 Li et al. (2007 , LCH) 
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Fig. 8. Synthetic spectra (continuous curve) and pure Pru line profiles (dashed curve, for /14223 A and /14449 A lines only) calcu- 
lated with the ATLAS9 solar model atmosphere in comparison with the observed NSO solar flux spectrum (Kurucz et al. 1984, bold 
dots). For the 4449 A blend, dotted lines show the synthetic spectra computed with +0.05 dex change in the adopted Pr abundance. 
The Prn lines are treated based on non-LTE line formation. 



-9.4 + 0.2 from 10 Pr n lines and by 2. 1 dex larger, log (Pr/H) = 
-7.3 + 0.3, from 14 Prm lines. Such a large discrepancy be- 
tween the rare-earth LTE abundances derived from two ioniza- 
tion stages is typically observed in the roAp and cool Ap stars 
(the REE anomaly, see Ryabchikova et al . 120041 1. 

Our non-LTE code does not include the magnetic field in 
statistical equilibrium calculations. The presence of a mag- 
netic field should not cause significant changes in the derived 
level populations, for moderate field strengths. As was shown 
in Sect.|U departures from LTE for the level populations are 
mainly caused by UV continua and by strong UV transitions. 
The change in emergent continuum flux around 2000A due to 
the influence of the magnetic field on atmospheric structure does 
not exceed 7%, for field strengths up to 10 kG (Kochukhov et al. 
2005). The critical strong UV transitions have small Zeeman 
splittings, since the splittings scale with the wavelength. Further 
abundance determinations are, therefore, performed in the field- 
free approximation. 



With fixed level populations, the magnetic field can affect 
an abundance analysis through Zeeman splitting of spectral 
lines. Therefore, we investigate a possibility to replace magnetic 
broadening effects by other broadening mechanisms, for in- 
stance, by a microturbulence. Different spectral lines used in this 
study have different Zeeman patterns. The Prn lines are weak 
and HFS-affected. We adopt a common value, V m ; c = lkms~', 
in their analysis. For each line of Pr m, a pseudo-microturbulence 
velocity is calculated as follows. We first fit the line profile with 
the SYNTHMAG code to derive the Pr abundance in LTE ap- 
proximation. Then, calculations are made with the SYNTH3 
code (Kochukhov 2007) that ignores the existence of magnetic 
field. The element abundance is fixed and V m i c varies until the 
theoretical equivalent width reaches the observed one. The indi- 
vidual values of V m j c are given in Table|4] Applying this proce- 
dure we find only minor changes in the element abundances, by 
0.14 dex and 0.16 dex from the Prn and Prm lines, respectively, 
compared to those determined from magnetic spectrum synthe- 



sis. It is worth noting, that the abundance difference between two 
ionization stages remains to be almost the same. 

Next, we check whether the departures from LTE can solve 
the problem of the praseodymium abundance in HD 24712 and 
find that non-LTE tends to reduce the difference in LTE abun- 
dances but fails to remove it completely. The non-LTE abun- 
dance is log(Pr/H) = -9.22 + 0.19 from the Prn lines, and a 
1 .90 dex larger value is determined from the lines of Pr m. Thus, 
we cannot obtain consistent abundances of the praseodymium 
from the lines of two ionization stages assuming the uniform 
element distribution in the atmosphere of HD 24712, either con- 
sidering non-LTE line formation, or taking into account the split- 
ting of the spectral lines in magnetic field. We suggest, therefore, 
that similarly to the neodymium (Mashonkina et al. 2005) the 
praseodymium is distributed non-uniformly in the atmosphere 
of HD 24712, although the properties of the Nd and Pr distribu- 
tions may be different. 

6.2. Vertical distribution of Pr in the atmosphere 

Using the simplified step-function approximation as an initial 
guess for the Pr abundance profile, we modify it by the trial- 
and-error method based on non-LTE line formation trying to 
fit the observed equivalent widths of both Prn and Prm lines 
in the spectrum of HD 24712. Due to certain limitations of the 
present non-LTE calculations, mainly due to ignoring Zeeman 
splitting, we cannot use the observed line profiles in stratifica- 
tion analysis as this was done, for instance, in LTE analysis for 
Ca (Ryabchikova et. al 2008). The final solution is shown in the 
top panel of Fig. [7] The praseodymium is concentrated mostly 
above log T5000 = -4, and the required Pr abundance in the layer 
is [Pr/H] > 4. TableH (column 10) presents the corresponding 
non-LTE equivalent widths of the Prn and Prm lines. For com- 
parison, we give also the LTE equivalent widths, calculated with 
the obtained Pr distribution, and the non-LTE abundance correc- 
tions. It can be seen that net non-LTE correction to the Pr abun- 
dance Anlte(Pi" ii) - ANLTE(Pr III) is at the level of 1 .8 dex, for the 
stratified atmosphere. In order to make the LTE element abun- 
dances from the lines of Prn and Prm consistent, the layer of 
enhanced praseodymium has to be located above the uppermost 
depth point in our model, log T5000 = -8 (!). This means that any 
stratification analysis for the praseodymium in roAp stars has to 
be performed based on non-LTE line formation. It is worth not- 
ing that a location of the enriched praseodymium layer turns out 
to be very similar to that found for the neodymium in our earlier 
study. 

A quality of fitting the observed praseodymium lines in 
HD24712 with the derived Pr stratification profile is illustrated 
in Fig. [9] and Fig.[l0] The synthetic flux profiles are convolved 
with the instrumental profile (spectral resolution R = 80000), 
projected rotational velocity v e sin / = 5.6 kms 1 and additional 
broadening parameter equivalent to radial-tangential profile of 
V ma c= 4kms~' for the Prn and 6kms~' for the Prm lines. One 
cannot expect to get the same quality fits for the magnetic star 
with a complex chemically stratified atmosphere as those for the 
Sun. Nevertheless, our non-LTE stratification analysis allows to 
reproduce reasonably well the observed Prn and Prm lines in 
HD 24712, contrary to the case of the uniform distribution of 
Pr. In the latter case, we can fit the Prm 5300A and 6195A 
profiles (dashed curves in Fig. [Tot with the Pr abundance deter- 
mined from non-LTE analysis of the Pr m lines only, log(Pr/H) 
= -7.32, but calculate extremely strong lines of Prn compared 
to the observations. 
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The knowledge of the Pr abundance distribution is impor- 
tant for a study of atmospheric pulsations. Fig. QT| shows a dis- 
tribution of the pulsation RV amplitudes and phases of the Pr 
lines as well as of the lines of Ndn and Ndm and of the Ha 
core in the atmosphere of HD 24712. The RV data are taken 
from Ryabchikova et al. d2007bl Online Table 4). We use the 
average optical depths of line formation which are calculated 
with the element abundance distributions obtained in this study 
for the Pr and by Mashonkina et al. (2005 ) for t he Nd . We fol- 
low the formalism suggested by Achmad et al. ( 1991 ) with us- 
ing the contribution function to the emergent line radiation. The 
computed values log t' 50QQ are given in Table [4] (column 12). For 
comparison, we present there (column 13) also the average op- 
tical depths of line formation based on the contribution function 
to the emergent total (line + continuum) radiation, logr^ . It 
was emphasized by Achmad et al. that the first approach pro- 
vides more realistic depth of line formation, in particular, for 
weak lines, for which the continuum contribution shifts an aver- 
age depth formation downward. This effect becomes much more 
important for the case of stratified atmosphere. When the con- 
tribution function to the total radiation is used, the formation 
depth of most Pr n lines is shifted well below the enhanced Pr 
abundance layer. The non-LTE formation depths for the Ha core 
have been computed according to Mashonkina et al. ( 2008J. It 
can be seen from Fig.Qj] that with the stratified Pr and Nd abun- 
dance distributions found based on non-LTE line formation, we 
get a consistent picture typical for running pulsation wave where 
the amplitudes and phases grow towards upper layers. A gap of 

~0.5 dex at log T5000 4.5 may be caused by the uncertainties 

of the modelling. 

6.3. The uncertainties of Pr stratification analysis 

In this subsection, we discuss the effects due to the uncertainties 
of atomic parameters and due to ignoring magnetic intensifica- 
tion of spectral lines. 

Photoionization cross-sections. The statistical equilibrium 
of Pr ii/m in the stratified atmosphere is mainly defined by en- 
hanced photoionization of Pr n. We perform test calculations us- 
ing various photoionization cross-sections. When the hydroge- 



Table 4. Equivalent widths (in mA) of the Pr n and Pr m lines observed in HD 24712 (column 4) and calculated using non-LTE and 
LTE approach for the model 7250/4.3/0 with the homogeneous (columns 6-8) and stratified (columns 9-11) element distributions. 
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nous photoionization cross-sections are increased by a factor of 
30 the overionization of Pr n is amplified and Anlte increases by 
0.03 dex to 0.05 dex for different Pru lines. The uncertainty of 
photoionization cross-sections is much less important for Pr m. 
Reducing the hydrogenous photoionization cross-sections by a 
factor of 100 has much larger effect on both the Pru and Prm 
lines. In this case, the layer of enhanced Pr has to be shifted out- 
ward by AlogT5ooo - 0.5 to agree the abundances of Pr from 
two ionization stages. The obtained Pr abundance distribution is 
shown by dashed curve in the top panel of Fig. [7] We empha- 
size that even with the lowest photoionization cross-sections the 
ionization equilibrium Pr n/Pr m deviates significantly from the 
TE one and non-LTE removes approximately 1 .5 dex of the dif- 
ference between the Pr abundances derived from the Pru and 
Pr in lines under the LTE assumption. As was discussed in our 
earlier paper (Mashonkina et al. 2005), the adopted hydrogenic 
approximation gives, probably, the low limit for the photoioniza- 
tion cross-sections for the levels in the rare-earth elements. So, 
we under- rather than overestimate non-LTE effects for Pr ii/m 
in our calculations. 

Collision rates. In SE computations, we assume an effec- 
tive collision strength to be equal T = 1 for every forbidden 
transition. This may be incorrect. The 7?-matrix method calcu- 
lations for electron impact excitation in Can (Melendez et al. 
2007) show that T depends on the transition energy separation. 
For example, at electron temperatures we are concerned with, 
T > 30 for the forbidden transitions with AEn < 0. 1 eV, while 



T < 0.1 if AEij > 10 eV. For the fine-structure transitions in 
Fen, the /^-matrix method predictions of Ramsbottom et al. 
(2007) give Y > 5. We perform non-LTE calculations for Prii/m 
assuming that for the forbidden transitions, Y depends on the 
transition energy separation as follows: log Y =1.5 for AEy < 
0.1 eV, log Y = -0.38 AE U + 1.54 for AE U = 0.1 eV - 4 eV, 
log Y = -0.17A£, 7 4- 0.67 for AE U = 4 eV - 10 eV, and log Y 
= -1 for AE^ ^ 10 eV. These approximations are based on an 
extensive set of the data of Melendez et al. (120071 1 for the for- 
bidden transitions in Ca n. With stronger collisional coupling the 
low-excitation levels to the Pr n ground state the departures from 
LTE are weakened for the Pru transitions arising from the levels 
with £e XC < 0.8 eV, and they are strengthened for the remaining 
transitions in Pru and for all transitions in Prm. The Pr abun- 
dance distribution obtained from these test calculations is shown 
by dotted curve in the top panel of Fig. [7] We find that no re- 
vision is required outside logrsooo = -6.5 and inside logrsooo 
= -4 compared to the distribution found in Sect. 16.21 while the 
praseodymium abundance is decreased by up to 0.2 dex in be- 
tween. 

Magnetic field. The magnetic field effects are taken into 
account in this study approximately by introducing a pseudo- 
micro turbulence. First, this works only for strong spectral lines. 
Second, magnetic field and microturbulence affect the saturated 
line profile in different ways. Magnetic field splits the line into 
the Zeeman components making the line less deep (magnetic 
desaturation) and broader, while microturbulence does not affect 
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Fig. 10. A comparison between the observed (bold dots) and non-LTE spectra for the selected Prn and Prin lines in HD 24712. The 
non-LTE profiles from calculations with the derived Pr distribution are shown by continuous curves. Dashed curves correspond to 
the non-LTE profiles from calculations with log(Pr/H) = -7.32 everywhere in the atmosphere. 



the line depth, but increases the half-width of the line profile. 
As was shown in Sect. 16. II taking into account Zeeman splitting 
produces small changes, by 0.16 dex, on average, in the derived 
Pr abundances compared to those determined without magnetic 
field. Such an effect cannot influence our conclusion on strati- 
fied distribution of Pr in the atmosphere of HD 24712. Zeeman 
splitting may also result in different line formation depths com- 
pared to that for turbulently broadened lines and, therefore, in a 
different location of the Pr enriched layer. We estimate this ef- 
fect using the strongest Prm /15300A line. Zeeman splitting of 
this line is represented by the three components. From compar- 
ison of the computed single line with V m \ c = 1.05 kms~' with 
the triplet with zero microturbulence we deduce that the mag- 
netic desaturation may lead to only a 0.2 dex decrease of the Pr 
abundance in the layer and in the shift of Pr abundance profile 
downward, by 0.35 dex in log T5000 scale. 

Thus, the uncertainty of the used atomic data on photoioniza- 
tion cross-sections and electron impact excitation cross-sections 
and ignoring the magnetic field in non-LTE analysis cannot 
change our conclusion as regards an existence of inhomoge- 
neous distribution of Pr in the atmosphere of HD 247 12. 

7. Conclusions 

In this study, we first presented a comprehensive model atom for 
Prii/in based on the measured and the calculated energy levels. 
Non-LTE line formation for an extended list of Prn and Prm 
lines was considered for the Sun and for the temperatures char- 
acteristic of A type stars, T eS = 7250 K - 9500 K. At r eff < 
8000 K, Prn represents the majority of the element and the de- 
partures from LTE for the Pr n lines are small. They are mainly 
caused by a deviation of the line source function from the Planck 



function and may result in the non-LTE abundance corrections of 
different sign for various lines. Overionization of Prn at T e g > 
8500 K leads to depleted total absorption in the Prn lines and 
positive Anlte that grows with temperature. Non-LTE leads to 
strengthening the Pr in lines and negative abundance corrections 
over the whole range of stellar parameters. As Prm becomes the 
majority species, Anlte of the Prm lines decreases in absolute 
value. 

Using the modern laboratory data for oscillator strengths of 
the Prn lines, improved partition function of Prn based on the 
laboratory and the calculated energy levels, we determine the 
non-LTE abundance of praseodymium in the solar atmosphere, 
log(Pr/H) G = -11.15 + 0.08 and log (Pr/H) = -1 1.10 + 0.08, 
for the ATLAS9 and the Holweger & Miiller (119741) models, re- 
spectively. The Pr abundance based on the theoretical model is 
0.05 dex to 0. 10 dex larger compared to the meteoritic value rec- 
ommended by Grevesse et al. (1996), Palme & Jones (2005 ), 
Lodders (2003 ), or Asplund et al. (120051) . but the deviation is 
still within the error bars. 

For the roAp star HD 24712, the element abundances from 
two ionization stages, Pr 11 and Pr m, reveal a discrepancy of two 
orders of magnitude if the atmosphere is assumed to be chemi- 
cally homogeneous. Neither non-LTE, nor magnetic effects can 
be responsible for such abundance anomaly. 

Introducing the layer with a strongly enhanced Pr abundance 
in the outer atmosphere provides a natural possibility to describe 
the lines of both ionization stages, Pr 11 and Pr m, for the sin- 
gle element distribution. We find that the required Pr overabun- 
dance in the layer is [Pr/H] > 4 at log T5000 < -4. The obtained 
praseodymium and neodymium stratifications provide a possi- 
bility to explain the distributions of the pulsational characteris- 
tics (radial velocity amplitudes and phases) over the significant 
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abundances derived from the lines of two ionization stages in the 
roAp stars. As was shown for Nd by Mashonkina et al. (2005) 
and for Pr in the present work, the departures from LTE for the 
lines of the first and the second ions are of the opposite sign, and 
they are large if the element is concentrated in the uppermost 
atmospheric layers where collisions are inefficient to establish 
thermodynamic equilibrium. In such a case, a stratification anal- 
ysis for REE has to be performed based on non-LTE line forma- 
tion. 
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Fig. 11. Pulsation radial velocity amplitude (top panel) and 
phase (bottom panel) variations for the lines of Pru (open di- 
amonds), Prm (filled diamonds), Ndn (open triangles), Ndm 
(filled triangles), and the Ha core (filled circles) as a function of 
the optical depth in the atmosphere of the roAp star HD 24712 
with the stratified distribution of Pr and Nd. 



part of the atmosphere of HD 24712, where the lines of these 
elements as well as the core of the Ha line are formed. 

Similar abundance profiles found empirically for the Nd 
(Mashonkina et al. 120051 1 and the Pr in the atmosphere of 
HD 247 12 point to some common physical mechanisms produc- 
ing accumulation of these two elements in the uppermost atmo- 
spheric layers. The main mechanism can be radiatively driven 
diffusion. The theoretical diffusion calculations of LeBlanc & 
Monin d2004t for Cr and Fe and of Alecian & Stift d2007i l for 
Mg, Si, Ca, Ti, Fe predict the existence of the abundance jumps 
up to several orders of magnitude in stable stellar atmospheres 
with effective temperatures characteristic of Ap stars. Similar 
abundance profiles are derived empirically for these elements in 
the atmospheres of roAp stars (Ryabchikova 2008). Opposite to 
REEs, Mg to Fe are accumulated in the lower atmosphere. Due 
to difficulties with the atomic data no diffusion calculations have 
been performed for the REEs to check the validity of the empir- 
ical distributions, but this situation is improving and we hope to 
see first REE diffusion calculations in the future. 

We are planning to extend non-LTE line formation study 
to other REE which reveal a discrepancy between the element 
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Fig. 1. The Pru even term structure. 
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^ I I I I II I I I I Mil II III 

-3 inn i mi i ii ii inn ii iii ii nil iiiiiiiiiiiiiiiiiiiiiii 

4 I III III I Ml II HIM llllllllllllllllll IIIIIIIIIIIIIIIIIIIIIIIIIIIIIH 

1 iii iiiii mi mi i ii ii mil i iiiiiiiiiiiiiiiiiiiiiii iiiiHi 

q 1 1 IIIIHI I III II IIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIIHH 

4 llll MM llllllllll lllllllllllllllllllllllllllllllllllllllllllimillll 

-o iiiii iii iii nun iiiiiiiiiiiiiiiiiiiiiii ii iii iiiiiiiiiiiiiiniimi 

* Ml III I IIIII llllll III II III I Ml I II IMIIIIIIIIIIIIIIIIIIIIIIIHII 

r-r in iiiii linn 1 1 mum nun i iiimi imiii limn 



"2 I II I II II II I I I I 

i in ii nil mi i iiimi iiii i mini i 

iiiii iii iiiii ii i iiiiiiiiii i iiii illinium i iiiiiiiiiiM 

1 mi iiiii iiiii iiiii illinium iiiiiiiiiiiiiiiiiiiiHim 
o iiiiiiiiii mi mi i iiiiiiiiiiiiiiiiii mi i illinium ii«h 
4 iiiii iiiii iiiii iniiiiiiiiiiiiii iiiii in mil iiiiiiiiiiiiiiiiii 

IIIIIIIIIIIIIIIIII IIIIIIIIIIIIIIIIIIIIIHIIII 

I IIIIIIIIII IIIII IIII IIIIIIIIIIIIIIIIIIIIIII 

IIII 



Fig. 2. The Pr n odd term structure. 
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Fig. 3. The Pr m model atom. 



